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We present detailed theory for Andreev level qubit, the system consisting of a highly transmissive 
^ ' quantum point contact embedded in a superconducting loop. The two-level Hamiltonian for Andreev 

levels interacting with quantum phase fluctuations is derived by using the path integral method. We 
also derive kinetic equation describing qubit decoherence due to interaction of the Andreev levels 
with acoustic phonons. The collision terms are non-linear due to fermionic nature of the Andreev 

CI 

states, leading to slow non-exponential relaxation and dephasing of the qubit at temperature smaller 
than the qubit level spacing. 
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I. INTRODUCTION 



The possibility to employ Andreev bound levels in superconducting contacts for quantum computation has been 
suggested in Refs. 1, 2. The proposed Andreev level qubit (ALQ) consists of a highly transmissive, with reflectivity 
R <C 1, quantum point contact (QPC) embedded in a low inductance superconducting loop. In the ALQ, quantum 
' information is stored in the microscopic two- level system of Andreev bound states in the contact. Hybridization 
of the clockwise and counterclockwise persistent current states in the ALQ loop is produced by the microscopic 
processes of electronic back scattering in the QPC. This is different from the macroscopic superconducting flux 
qubits 3-5 and charge-phase qubit 6 , where the hybridization is provided by charge fluctuations on the tunnel junction 
capacitors. Thus requirement of the large charging energy or large loop inductance is not critical for the ALQ. A 
fNJ \ single ALQ consists of a pair of the Andreev bound levels belonging to the same normal conducting mode in the QPC; 

a multimode QPC will form a qubit cluster. The way of the ALQ operation is similar to the one of the experimentally 
| tested flux qubits 3-5 - the Andreev levels can be excited by driving biasing magnetic flux through the qubit loop. 7,8 ' 1 
The read out method is also similar to the flux and charge-phase qubits: the quantum state of the Andreev levels 
determines the magnitude and direction of the persistent current circulating in the loop, and also the magnitude of 
the induced flux. Since the quantum information is stored in the microscopic system, Andreev levels, while the access 
for manipulation and readout is provided by macroscopic persistent currents, ALQ occupies an intermediate place 
between the microscopic solid state qubits (like localized spins on impurities 9 or quantum dots 10 ) and macroscopic 
■ superconducting qubits. 

During 90-s, the Josephson transport in superconducting QPCs has been intensively investigated, and a num- 
ber of remarkable experiments has been performed on atomic size metallic QPCs using controllable break junction 
technique 11 ' 12 , as well as on gated quantum constrictions in 2D electron gas confined between superconductors. 13 The 
critical Josephson current and current-voltage characteristics have been thoroughly examined in these experiments 
by applying current or voltage bias. 13-16 There was, however, one experiment, particularly important in the qubit 
q . context, where the flux bias has been implemented: Koops et al. 17 have inserted metallic QPC in a SQUID, and eval- 
uated the Josephson current-phase dependence by measuring the induced flux with an inductively coupled SQUID 
magnetometer. The measurements have been only reported for the equilibrium state. Unfortunately, no experimental 
attempts to drive the QPC out of equilibrium to some coherent or incoherent excited state have been performed so 
far. 

The purpose of this paper is to present a theory for the Andreev level qubit in more detail than that was possible 
in short communication. 2 We will also consider the electron-phonon interaction as a potential "intrinsic" source of 
the qubit decoherence, and derive a kinetic equation for the qubit density matrix. The ALQ Hamiltonian and the 
kinetic equation are derived by using a path integral method. 18-20 The central technical difficulty here is to extend the 
method to contacts with large transparency. This difficulty is overcome by incorporating exact boundary condition 
in the QPC action. Another important point discussed in the paper is the role of the charge electro-neutrality in the 
junction electrodes, which affects the qubit Hamiltonian. 

The fermionic nature of the Andreev levels does not affect the qubit operation and qubit-qubit coupling but it plays 
important role for the qubit decoherence. We find that the electron-phonon collision terms in the kinetic equation for 
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the ALQ differ qualitatively from the Bloch-Redficld master equation 21 commonly applied to study decoherence of the 
macroscopic superconducting qubits. 22 This results in long phonon-induced decoherence time for the ALQ, namely, 
at the temperature smaller than the qubit level spacing, both the relaxation and dephasing processes are governed by 
a power rather than exponential law. 

The structure of the paper is the following. We discuss the model description of the QPC in Section II, and then 
explain in detail, in Section III, the path integral approach for transmissivc QPC: we consider the action for the 
contact and derive the effective Hamiltonian for Andreev levels interacting with the quantum phase fluctuations; the 
single-particle density matrix, and effective current operator are also discussed in this section. In Section IV we discuss 
averaging over fast phase fluctuations and derive an effective Hamiltonian for the qubit, this procedure is extended 
in Section V to the two inductively coupled qubits to derive direct qubit-qubit interaction. Section VI is devoted 
to the electron-phonon interaction: we derive effective action for the Andreev level-phonon interaction and derive 
corresponding collision terms in the equation for the qubit density matrix; then we present solution of the kinetic 
equation, and evaluate the decoherence rate. 

II. CONTACT HAMILTONIAN 

Let us consider a superconducting quantum point contact with bulk 3D electrodes. We model the contact with a 
smooth on the atomic scale (adiabatic) constriction and assume a local scatterer situated in the neck of constriction 
(see Fig. 1) producing weak electronic back scattering with the reflectivity iJ« 1. We further assume that the 
constriction supports a single conducting mode. 
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FIG. 1. Adiabatic superconducting constriction with a local scatterer (dark region) in the neck. The length of the constriction 
is small on the scale of the superconducting coherence length but large on the atomic scale. 

We adopt the mean field approximation for the electrons in the contact, which are described with the Hamiltonian, 

H e = j dv *t( r , t )h(r, t)*(r, t) + \ CV 2 {t), (1) 

where the first term is the BCS Hamiltonian for the bulk superconducting electrons, ^(r, t) being the two-component 
Nambu field operator, and the second term describes the charging energy of the contact capacitor C 18 20 . The 
single-particle Hamiltonian h in Eq. (1) has the form, 

a z + A(r,t)e^^a x , (2) 

where A(r,i) and x(r,i) are, respectively, the modulus and phase of the superconducting order parameter, the 
potential U(r) accounts for the confinement of electrons within the contact as well as the electron scattering, while 
<fi(r, t) and A(r, t) are electromagnetic potentials. The voltage drop at the contact, V(t), in Eq. (1) is related to the 
discontinuity of the electric potential at the contact, V(t) = (p(— 0, t) — ^>(+0, t). To investigate the decoherence effects, 
we allow the electrons in the electrodes to interact with acoustic phonons, and include corresponding electron-phonon 
interaction and phonon terms in the total Hamiltonian of the contact, 

H c = H e + H e - p h + H. p h ■ (3) 

It is convenient to introduce quasiclassical (Andreev) approximation for the superconducting electrons. Following 
standard procedure, we eliminate rapidly varying in space potential, U (r), by introducing quasiclassical wave functions 
of the single conducting mode in the left and right electrodes, 



h = 



HfiV-(e/c)A(r,f) a z y 
2m 



- fi + U{r) + eip{r, t) 
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tt(r,t) = ]T ip±(r±,x)e ia f dxk{x) ^\x,t), (4) 



cr=± 



and couple these wave functions by means of a normal electron scattering matrix. In Eq. (4), ip("\x) are slowly 
varying ID envelopes for the longitudinal electron motion (<r = ± indicates the direction of the motion) , ipj_ (rj_ , x) is 
a normalized wave function of the transverse motion with energy E±(x), p{x) — Kk(x) — \J2rn (jj, — E±(x)) — mv(x) 
is the quasiclassical longitudinal electronic momentum. The coupling of the quasiclassical envelopes in the left (L) 
and right (R) electrodes, tpL,R, is conveniently described by the transfer matrix, 7 ' 23 



Here d and r are the energy independent transmission and reflection amplitudes, respectively. Since any observable 
quantity is expressed through a bilinear combination of the envelopes with the same a, the energy independent 
scattering phases can be eliminated from the boundary condition (5), hence, without loss of generality, the scattering 
amplitudes will be further assumed to be real, d = \/D, r — \/R, where D and R are the transmission and reflection 
coefficients of the contact, respectively. 

The electromagnetic potentials, tp(r, t), A(r, t), and the complex order parameter in Eq. (2) are to be found from the 
Maxwell equations and the self-consistency equation. It is convenient to present the Hamiltonian in a gauge invariant 
form by extracting the phase of the order parameter using a local gauge transformation, ^(r, t) — ► e lCT *x(r,*)/2 ^ 
Then the superfluid momentum p s = ?iVx/2 — (e/c)A, and the gauge-invariant electric potential (p = h\/2 + eip, 
appear in the quasiclassical Hamiltonian of the electrode, 

= av{-i%d x ) a z + avp s + (pa z + Acr x , (7) 

while the phase difference <j>(t) = Xfi(0,i) — Xi(0,t), appears in the boundary condition, 

J> _t e icr z 4>(t)/2 ji _ ^ 

In the bulk metallic electrodes with good screening, and at small frequencies relevant for the problem, the 
gauge-invariant fields, tp(r,t), p s (r, t), are to be found from the electro-neutrality condition, and the current 
conservation, 24-26 

5n(T,t)=0, Vj(r,t) = 0, (9) 

where n(r, t) is the electronic density. In the electrodes, the charge imbalance relaxation yields the equilibrium relation, 
ip = (dn/dfj,)~ 1 Sn, on the distance exceeding the electric field penetration length. Furthermore, in the absence of 
normal dissipative current, p s is proportional to the total current density, j = (en/m)p s , which is negligibly small far 
from the contact due to rapid spreading out of the current in the point contact geometry. Thus the conditions (9) 
yield complete cancellation of the electromagnetic potentials in the electrodes, 

<p(x,t) = 0, p,(x,t)=0. (10) 

Taking into account that the modulus of the order parameter far from the contact is equal to the equilibrium value, 27 
A = const, we conclude that the Hamiltonian (7) approaches the equilibrium form. 

Potential (p can be expressed through the electric field E, (p = hx/2 + e J dx E, by introducing the gauge invariant 
phase x — X ~ (2e/cK) J dx A — (2/h) J dxp s . The spatial distribution of these quantities is illustrated on Fig. 
2. Since (p vanishes far from the contact, and \ is constant according to Eq. (10), then the gauge-invariant phase 

difference across the contact 4> = Xh(0> — Xl(0> is related to the voltage drop, V = dx E , 

kt) = ^ (ii) 

(Josephson relation). In the SQUID, this relation is equivalent to the phase versus flux relation, <j> = 2e$/fic, since the 
voltage drop across the contact is generated by time variation of the magnetic flux $ threading the SQUID. Finally, 
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we notice that the gauge- invariant phase difference, <f>, rather than <f> enters the boundary condition (8), which can 
be explicitly seen by extracting the Aharonov-Bohm phase from the transfer matrix. Thus the gauge-invariant phase 
difference remains the only free collective variable whose dynamics is determined by the electrodynamic environment 
of the contact. Below we will not distinguish between 4> and <j>, because the difference is negligibly small in the QPC. 




FIG. 2. Spatial distribution of the electric potential, ip (thin line), gauge invariant potential, <p (bold line), and the time 
derivative of the gauge invariant phase, \ (dashed line), in the vicinity of the contact. 



Proceeding with discussion of the interaction of electrons with phonons, we consider only longitudinal acoustic 
phonons and describe the interaction within the deformation potential approximation, 

H e - ph = 1 j drVu(r,t)* t (r,t)a z *(r,t), (12) 

where 7 is the deformation potential constant, u(r, t) is the phonon field operator, 



»w = E - q ^ eiqr + b i e ^ qr ) ■ = ( 13 ) 

s is the sound velocity, and p is the crystal mass density. The Hamiltonian of free phonons has a standard form, 

H ph = m « ( b l b <* + V2) • (14) 
q 

Our strategy now will be to derive effective Hamiltonian for the Andreev levels including interaction with phonons. 
If the phase difference would be a classical variable, this derivation can be in principle done by direct truncation of the 
Hamiltonian (3). However, in the presence of quantum phase fluctuations it is convenient to apply the path integral 
technique. 



III. CONTACT EFFECTIVE ACTION 



Let us consider the whole system, the QPC and superconducting loop (see Fig. 3), and introduce the path integral 
representation for the propagator, 

U = jv 2 iP L V 2 i; R V{Xc l }V(be l I dtLt ° t/n , X q = 6 q + &* . (15) 

The Lagrangian of the system L to t consists of the contact part L Cl and the part describing the circulating current in 
the loop. The latter is conveniently combined with the charge term from the electronic Hamiltonian (1) giving the 
Lagrangian of the loop oscillator L OSCl 

L tot = L C + i osc , L osc - (J£) (d t( f>) 2 " ^ (0 " 0e) 2 ) , (16) 

where <j) e corresponds to the bias magnetic flux, and L is the loop inductance. The remaining part of the contact 
Lagrangian consists, similar to Eq. (3), of the electronic part, phonon part, and electron-phonon interaction, 

L c = L e + Lph + Le—ph . (17) 
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In the quasiclassical approximation, the electronic Lagrangian splits into two parts, L a , a = L,R, corresponding 
to the left and right electrodes, 

L a =J2 I dx t>* } ( x ' *) ( x ' *) ^ ( x > *) ' ( x ' *) = iW t + ™iKd x a z - Aa x , (18) 
<j=± •> 

and the third part Lbc discussed in detailed in the next section, which accounts for the boundary condition. Noting 
that relaxation processes are caused by phonons with small wave vectors compared to the Fermi wave vector, q <C k, 
transitions between the states and are forbidden, and the electron-phonon Lagrangian can be written on 
the form, 

L e - ph = - 1 ]T Yl f dr H±\ 2 ^ a a) (x,t)a z ^Hx,t)Vu(r,t). (19) 

a=L,Ra=± J 

The phonon Lagrangian is given by equation, 

q 




FIG. 3. Sketch of the Andreev level qubit: a low inductance superconducting loop with a quantum point contact (QPC). $ 
is the magnetic flux; the arrows indicate fluctuating persistent currents. 



A. Boundary condition 

The boundary condition (5) is valid for any contact transparency. To include this boundary condition in the path 
integral formulation, we introduce an additional term in the Lagrangian 2 , 

Lbc = V(t) [^ e^'W^? (0, t) - (1 + VRy*'*™^ (0, t)] + h.c. , (21) 

where 77 is an auxiliary fermionic Nambu field playing the role of a Lagrange multiplier. Correspondingly, the inte- 
gration over 77 is to be included in the propagator in Eq. (15), giving the following form for the electronic part of the 
propagator, 

U e = J V 2 V V 2 ^ L V 2 ^ R e l I 'dt(L L +L n +L BC )/h (22) 

Let us prove that such a Lagrangian indeed generates the boundary condition (5), (8). To this end, it is convenient 
to make a step back and restore a non-quasiclassical form for the fermionic field, 

^(x,t)=J2e^I kdx ^\x,t), (23) 

(7 

in the bulk part of the Lagrangian, 

L L + L R =j dxi>(x,t)C(x,t)ip(x,t), C(x,t) = ihd t - [(-h 2 /2m)d 2 + E ± - fi\a z - A a x . (24) 

<V — oo 
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The dynamic equations and the boundary condition result from the zero variation of the action S e = J dt(LL + Lr + 
Lbc) with respect to ipL,R and fj, 



or in the explicit form, 



SS e 8S e . . 

£(x,t)i> L (x,t) = -VdS(x) e^'Wriit) , 

C(x, t) 1> a (x, *) = (! + Vi?) <5(z) e --^W/ 4 v (t) , (26) 



and 



\/^e- <<T ^ ( * )/4 VL(0,t) = (l + Vi?)e^ ( ' )/4 ^(0,t). (27) 
Integrating Eqs. (26) over x in a small vicinity of x — yields the relations, 

(7i 2 /2m)<7^ x VL(0,i) = y/De^'+WSit) , 

(h 2 /2m) a z d x i> R {0, t) = (1 + \/#) e- iff '*W/ 4 »?(*) . (28) 
Then introducing again the quasiclassical envelopes and combining Eqs. (27), (28) with the quasiclassical relation, 

0*Vl,*(O, t)=ikJ2a 4^(0, t) , (29) 

<T = ± 

we get the boundary condition equivalent to Eqs. (5), (8). 

B. Effective action for Andreev levels 

Now we are prepared to derive an effective action for the Andreev levels. Following the procedure of Rcf. 18, we 
integrate out fast electronic fields tp a in Eq. (22), 



(30) 



a—L.R<7—± ^ 

Then we are left with the effective action S°, which contains only variables, n{t) and 4>{t), 

S l = - J dt i dt 2 V(ti) [De-^^/V*! ~ *2) e la ^ (t2)/i + (1 + Vfi) 2 e la ^ {tl)/A g(ti - t 2 ) e -^-*(*0/4] v ( h ) , (31) 
where g(t) is given by the Fourier component, 

= fiUJ + A(T X 

9uJ hv^A 2 - (hw) 2 ' 

Connection between the effective action (31) and the Andreev levels can be established by considering the case of 
time independent phase, <fi = const. Indeed, by writing the effective action in the Fourier representation, 



_ 2(1 + VR) f duj 



/ q> , — q> 
Tiuj + A ( cos — a x — v R sin — a y 



Vu, ■ (33) 



we identify the spectrum of the system, huj = ±E a ((f>), by calculating eigenvalues of the matrix in the brackets, 

E a {4>) = A v /cos 2 (0/2)+i?sin 2 (0/2); (34) 

this equation coincides with well known Andreev level spectrum 28 ' 29 (see Fig. 4). Thus we conclude that the fermionic 
field n represents the Andreev levels. 
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Proceeding to a time-dependent 4>{t), we restrict its time variation speed to small values, Tidt4>/A <C A. Furthermore, 
the dynamics of the Andreev levels is also to be slow, E{(f>) <C A, which implies that the contact reflectivity must be 
small, R <C 1, ALQ must be biased at <p e w n, and the amplitude of the quantum phase fluctuations <fi(t) = <j>{t) — <f> e , 
must be sufficiently small, <f> -C <fi e - We emphasize that the constraint on the amplitude of the phase fluctuations 
is actually provided in our case due to the loop geometry of the electrodes having sufficiently small inductance; 
the constraint is important to suppress the Landau-Zener transitions to the continuum states. Under the imposed 
conditions, the non-local in time kernel g(ui), Eq. (32), can be replaced by a constant value (adiabatic approximation), 
^/A 2 — (fiuu) 2 — ► ^/A 2 — E% = ( e with E a = E{<j> e ), leading to the equation, 

2(1 + Vr) 



J dt r](t) ifidt + h<p a z + A ^cos ^ a x — Vr sin ^ a^j 



V(t) ■ (35) 



We eliminate the term with the phase time derivative in Eq. (35) by transforming 77, 



,2(1 + y/R) 

Then we finally arrive at the effective action 



^e_\ 12 e ia z VR4>/^ e i^/4 v _ (36) 



5° = J dtfj(t)[ihd t -h a }r)(t), (37) 

where 

r= . / rh _ rh \ 

(38) 



h a = A e --- V5*/2 ^ cos I Gz + ^ sin I ^ 



presents the effective single-particle Hamiltonian for the two-level Andreev system. 2 The Hamiltonian in Eq. (38) 
differs by the exponential pre-factor from the two-level Hamiltonian derived in Refs. 30, 31 and further employed 
in Refs. 32, 33. This factor appears in the present derivation after the electric potential has been included into 
consideration in Eq. (1) to provide the electro-neutrality in the electrodes [see text after Eq. (8)]. Both Hamiltonians 
are equivalent under stationary conditions, d t (f> — 0, and the difference is not important for the adiabatic dynamics. 
However, in general, the pre-factor is important, e.g. for derivation of correct equation for the current operator in 
Eq. (50). 




FIG. 4. Spectrum of the Andreev levels in a QPC with finite reflectivity (R — 0.04) (solid line), the level anticrossing is 
produced by electronic back scattering; at 7? = the Andreev levels (dashed line) coincide with the current eigen states. 

It is instructive to compare the case of the transparent contact considered here with the case of a tunnel contact 
extensively studied in the MQC theory 18-20 Physical difference between the two cases is that the Andreev level system 
in transparent contacts is a slow one, while in tunnel contacts it is a fast one because the Andreev level energy in 
tunnel contacts is close to A. Within the present formalism, the integration over ip fields is similar in both cases. 
However, the next step, the adiabatic approximation in Eqs. (31), (32) is not allowed in the tunnel limit; instead 
one should perform also the integration over r\ in Eq. (22), and make an expansion over small D. The result of this 
calculation, presented in Appendix A, coincides with the results of Refs. 18, 19, 34. 
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C. Andreev levels density matrix 



The Hamiltonian (38) governs the free evolution of the Andreev variable 77, 



SS° 

7^=0, ihd t r] = h a r]. (39) 



5fj(t) 



The same Hamiltonian governs the free evolution of the single-particle density matrix of the Andreev levels p a (t) 
defined via statistical average, 

Pa (t) = (mtf(t)), (40) 

where r) denotes the fermionic operator corresponding to the Grassmann field r], and the average is taken over all 
electronic states. The density matrix p a is a 2 x 2 matrix in the Nambu space, which satisfies the normalization 
condition Trp a = 1. The statistical average in Eq. (40) is represented, after the averaging over ^l,r, by a path 
integral, 

p a (t) = Jv^n{t)fi(t)e</ h . (41) 

Then calculating the time derivative, using Eq. (39) and the conjugated equation, we get, 

ihdtp a = [h a , Pa]. (42) 

Thus the Andreev level dynamics is described by a Liouville equation similar to ordinary quantum mechanical two-level 
systems. 



D. Andreev level current 



We conclude this section with derivation of an effective current operator for the Andreev levels. A common quantum 
mechanical expression for the current applied to the fields ipL,R at x — =f0 gives, 

h,R(t) - ^ [^.*(°' *) H^ x ) 1>l,r(P> *) + h - c ] > ( 43 ) 

or using Eq. (28), 

/(*) = IlA*) = (ie/h)VD fj{t) <j z e- l °^'^ L {Q, t) + h. c. (44) 

[the current is continuous, Ih{t) = lR(t), by virtue of Eq. (27)]. The same equation can be obtained by varying the 
electronic part of the action, S e , with respect to the phase difference, 

'M-™- < 45 » 

Using this relation the statistically averaged current, 

(/(*)) = J V\V 2 ^ L , R I{ty s *l n , (46) 

can be written as 

(J(t)) - 2e lj ^ f) J V 2 V V^ R e^' h , (47) 

or, after tracing out the fields, ipL,R, 

{I(t))=2ei-^jv\e iS VK (48) 
This leads to the following equation for the current in terms of the Andreev variable, 
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in which the kernel / corresponds to the effective single particle current operator of Andreev levels, 

'-%%—$H*)e->~<™> l ,.. K*)-AD*.f ,50) 

Apparently the current operator does not commute with the Andreev level Hamiltonian, [h a , I] 7^ 0, which is the 
consequence of a normal electron reflection at the QPC. Hence the Andreev level eigen states consist of superpositions 
of the current eigen states, unless R = (see Fig. 4). Correspondingly, the Andreev level current, defined as an 
expectation value of the current operator over the Andreev state, 

,2e dEgW) _ eDA 2 . 

/a = (/>a = ± ¥^^ T ^: sin ^ (51) 

differs from the current eigenvalues, ^feT/h [(...) a denotes averaging over the Andreev level eigenstate]. Thus the 
Andreev level current undergoes strong quantum fluctuations. The spectral density of current-current correlation 
function can be directly calculated by using Eqs. (38) and (50), [cf. Ref. 35], 

S(uj) = (I) 2 a Rt&n 2 ((f>/2) S(uj - 2E a /H). (52) 



IV. AVERAGING OVER PHASE FLUCTUATIONS 



Equation (42) describes dynamics of the Andreev levels for a given realization of the time dependent phase across 
the QPC. However, the phase dynamics is strongly coupled to the Andreev levels. Intrinsic dynamics of the phase is 
governed by the quantum Hamiltonian of the loop oscillator [cf. Eq. (16)], 

*™ = sb + -^. M] = ~ih, M=_, (53) 



where u> p = y SE^Ec/fi is a plasma frequency of the oscillator determined by the contact charging energy, Eq = 

e 2 /2C, and the loop inductive energy, El = (hc/2e) 2 (l/L). Thus the whole system is generally a multilevel one. On 
the other hand, the phase dynamics and the coupling between the loop oscillator and the Andreev levels is a vital 
part of the ALQ operation: The Andreev levels cannot be manipulated if the phase dynamics was frozen because any 
manipulation requires variation of the current, and second, the read out of the Andreev levels can be only performed 
via measuring the quantum state of the loop oscillator. An obvious way to solve this problem and preserve the qubit 
property of the coupled system is to "enslave" the loop oscillator by choosing the oscillator level spacing Hu> p much 
larger than the Andreev level spacing, tuv p 3> 2E a . Then the Andreev state evolution will not excite the oscillator, 
which will remain in the ground state and adiabatically follow the evolution of the Andreev levels. This implies that 
the phase is a fast variable which should be averaged out, leading to an effective qubit Hamiltonian. 




FIG. 5. Potential energy diagram for two displaced oscillator states corresponding to the different current states in the point 
contact (shown as arrowed circles) . The plasma frequency, lu p , is large compared to the Andreev level spacing; the oscillator 
remains in the ground state during the qubit evolution (dashed arrow) . 
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To facilitate the averaging procedure we take advantage of the small amplitude of the phase fluctuations, <p — 
4>— 4> e <C <j> e , which was already assumed when proceeding from Eq. (33) to Eq. (35). This assumption is justified by 
the large inductive energy, El 3> Ej w A, and it allows us to expand the Andreev level Hamiltonian (38) over small 
</>; then proceeding to the current eigen basis, r\ — -> e~ %<r " s ^' e l 4 "q, we g e ^ 

h a = A ^cos y g z + Vfisin y a y ^j - 0cr z = h° a + h int . (54) 

Averaging over fast phase fluctuations can be done directly in Eq. (15) by performing explicit integration over 
phase. However, there is a simpler way to get the same result. Equation (42) holds for a fluctuating phase provided 
the oscillator is not excited. It can be viewed as the diagonal with respect to <j> part of a more general equation for a 
full density matrix, p(a(j), a'cfi'), for the Andreev levels plus oscillator system, 

ihdtp = [h° a + h in t + H osc , p]. (55) 

The interaction between the Andreev levels and the oscillator given by the second term in Eq. (54) displaces the 
oscillator steady state from the origin by ±I/2M^ depending on the direction of the current in the junction or, 
equivalently, the state of the Andreev levels (see Fig. 5). We eliminate this term in (54) by applying the transformation, 

and then average the resulting Hamiltonian over the oscillator ground state, taking into account the relation, (e lAp ) = 
exp (— A 2 (p 2 } /2) [(...) indicates averaging over the oscillator ground state]. As the result we get an effective 
Hamiltonian, 



(h a ) Q = A |^cos — a z + VR* sin — a y J = h q (<t> e ), (57) 

where the bare contact reflectivity, R, is renormalizcd by the phase fluctuations, 

R* = e- 2X R, A = I 2 (c/) e )/4MhLu*. (58) 

This renormalization effect can be understood as the effect of inertia of the loop oscillator, which hinders the cur- 
rent variations, i.e. it works against the effect of the electronic back scattering at the contact responsible for the 
hybridization of the current states [see discussion in the end of the Section HID] . The renormalization effect becomes 
increasingly strong in the limit of a classical oscillator with large " mass" . Because of renormalization of the contact 
reflectivity, the Andreev level spectrum is modified, 

E a (^ e ) -> E* a (<j> e ) = A\J cos 2 (0 e /2) + R* sin 2 (0 e /2) , (59) 

and the qubit frequency reduces. This might be important for practical applications, because it would allow one to 
reduce the qubit frequency by choosing the circuit parameters rather than by tuning the contact reflectivity. 

Eq. (57) gives effective Hamiltonian h q for the ALQ in the absence of interaction with phonons. Similarly, averaged 
over phase fluctuation density matrix (p)o gives density matrix for the ALQ. Keeping the same notation, p, for the 
qubit density matrix, we finally arrive at the equation for the free evolution of the ALQ, 

ihdtp = [hq(<f> e ) , p]- (60) 

This equation is sufficient to describe the ALQ manipulation (by driving the biasing flux) and read out (by measuring 
the induced current) 1,2 , and also the qubit-qubit interaction, which is discussed in the next section. 



V. INDUCTIVE QUBIT-QUBIT COUPLING 



Our treatment of the interaction of the Andreev levels with phase fluctuations can be easily extended on the case 
of several inductively coupled SQUIDs to describe direct qubit-qubit coupling. Let us consider as an example two 
SQUIDs with different QPC reflectivities, R\ ^ R2, and identical circuit parameters, C and L, mutual inductance is 
A4. Then the inductance terms in the Lagrangian (16) written for the two qubits will take the form, 
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Us) ^ 



02 



i M 
M L 



(61) 



where 0i. 2 ar e fluctuating phases in the first and second qubits. By introducing the normal modes for the LC- 
oscillators, (j> — > e~ lTy7T / 4 (f), we present the two-qubit Hamiltonian on the form similar to Eqs. (54), (56), 



H = E 



i=l,2 



A cos — ^— <j z i + v Ri sin — ^ ( 



2\/2 



((-1) 4+1 X 1( T Z1 +X 2 a z2 ) + F 



(62) 



where H osc ^ describes the normal oscillator with frequency u> P i = cj \[C (L ± ~M), and the index i = 1,2 refers to 
the first and second qubit. Then we apply the similar transformations as in the previous section to Eq. (62), namely, 
we eliminate the linear in <j>i terms by means of a canonical transformation, H — > exp (i J2i Mvi) H exp {—i ^ AiPi), 
with Ai = ((— iy +1 lia z i + T 2 a z2 ) I (2^/2Mhujf), and then average the transformed Hamiltonian over the ground state 
of the normal oscillators. As the result, we obtain an effective two-qubit Hamiltonian including direct qubit-qubit 
coupling, 

H = h ql + h q2 + {e/hcf MZ{I 2 a 7A o z2 . (63) 

The renormalized contact reflectivities are now given by equation, R* = Ri exp[— (lf/4hM) (wf 3 + uj^ 3 )]. 

It is worth mentioning that the two-qubit configuration may be also realized with a single SQUID containing QPC 
with two conducting modes. In this case, we have the two Andreev level Hamiltonians coupled to the same loop 
oscillator, 



H 



E 



=1,2 



T i (0 e ) I 
h a,i -Z—9°z 



(64) 



Averaging over the phase fluctuations, we arrive at the same interaction Hamiltonian as in Eq. (63) but with the loop 
inductance L substituting for — M.. Thus we come to an interesting conclusion that the effect of the phase fluctuations 
not only reduces the bare contact reflectivity but also introduces effective coupling of the Andreev levels of different 
conducting modes in multimodc QPC. 



VI. ANDREEV LEVEL - PHONON INTERACTION 



A. Effective action 



In this Section we include electron-phonon interaction in the consideration. We start with the derivation of an 
effective action for the Andreev level-phonon interaction. To this end we repeat the calculation of the previous section 
adding the Lagrangian L e - p h, Eq. (19), to Eq. (30). By keeping in the electrode Green functions only the first-order 
correction in small interaction constant 7, we arrive at the following action, 

S v = - J dttdttfjih) [De-^^/ i G L {t 1 ,t 2 )e^ ( - t2 '>/ i + (1 + s/Rf e^^^GR^, t^e-™*^'^ V (t 2 ) , (65) 

where the Green functions Gl,r read, 



G a (t u t 2 )= 5i CT) (0,0;ti-i 2 )+ 7 / dr\^\ 2 f dt g£\0, x; h - t) Vu(r, t) a z g^ {x, 0; t - t 2 ) 



(66) 



In this equation, the quantities g^ R refer to the different parts, g^(x < 0,x' < 0;t) and g^(x > 0,x' > 0;t), 
respectively, of the translation-invariant free electron Green function, g^(x — x',t), which obeys the equation 



[ihd t + aa z vihd x — Aa x ] g^ a \x — x',t) = S(x — x 1 ) 8{t) . 



(67) 



The Green functions g^ in Eq. (66) are explicitly given by 
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gg,(x,0;t) = -Q(tx) 6 {ihdt + Aa x T<Ji(t<J z }6(t), (68) 

where Q is given by the Fourier component, = \/A 2 — {fiuo) 2 . 

Proceeding to the adiabatic approximation discussed in the previous section (£ w — > £ e ), an d performing the trans- 
formation (36), we arrive at the following effective action for the Andrccv lcvcl-phonon interaction, 



S v - P h = -J J dt J dr \tp±\ 2 n(x,t) Vu(r,i) , 



(69) 



where 



e -2\x\t e /hv 

n(x,t)= - — (1 + sgna; VR) A(x, t) , (70) 

and 

A(x, t) = fi 2 {d t fj) a x (d tV ) + TiA (d t f) e -to-(Vfl-gnx)0./2 o . w)J + ft c ^ _ S 2 - ^ 

Taking into account the zero-order dynamic equation with respect to 7, Eq. (39), and putting <f> = <j> e , 

i%d t r, = A e-f.SR*./* ^ cos ^ a z + VR sin ^ r? , (72) 

we obtain for the quantity n(x,t), Eq. (70), the following expression, 

n(x, i) = (-/R/2) k sgnx e^ 2 * 1 * 1 fja x ri, (73) 
where k = ( e /hv. Finally, integrating over r in the action (69), we present the effective interaction on the form, 

S v -ph = - / dt ^ 7q^q»? ^xV , (74) 

J q 

where 7 q is the effective constant of the Andreev level-phonon coupling, 



7q = 7«y^§^ " dxF( q± ,x)e- 2 ™sm(q x x), F(a ± ,x) = J dr ± |^(r)| 2 e^ r - . (75) 

We notice that the Andreev level-phonon interaction in Eq. (74) has purely transverse origin, i.e. while inducing 
interlevel transitions and hence the relaxation it does not produce any additional dephasing to the one associated 
with the relaxation. 

It is important to mention that the effective coupling constant is proportional to V^R, and it turns to zero in the case 
of perfectly transparent constriction. This results from the already mentioned fact that the relevant phonons have 
small wave vectors and are not able to provide large momentum transfer (~ 2hk) during scattering with electrons. 

The effective action in Eqs. (74) and (75) was derived for a given realization of the time dependent phase. To take 
into account the effect of phase fluctuation, we have to apply the transformation, Eq. (56), to the action (74); the 
Lagrange form of the transformation reads, 

ry-oxp^^^ fdrfafjr,. (76) 

Then the integration over the phase adds the factor e~ A to the action, cf. Eq. (58), which simply implies a renormal- 
ization of the coupling constant, 

7q - 7* = e- A 7q ■ (77) 

This result can be expected: since the Andreev level-phonon coupling is transverse in the current basis and depends on 
the contact reflectivity, the renormalized reflectivity R* rather than the bare one has to enter the coupling constant, 
7* = -f q (R*); this is consistent with Eqs. (77), (75). 
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B. Kinetic equation 



Qubit decoherence is usually described with collision terms in the Liouville equation for the qubit density matrix, 
which take into account interaction with an environment. Our goal now will be to derive the phonon- induced collision 
terms in the equation (60) for ALQ density matrix, and to evaluate the decoherence of the ALQ. 

While the description of free qubit evolution was possible in the terms of the single-particle density matrix, evaluation 
of the collision terms goes beyond the single-particle approximation and requires the knowledge of electronic two- 
particle correlation functions. This is because the Andreev levels are not a rigorously isolated system but rather belong 
to a large fermionic system of the superconducting electrons in the contact electrodes. Thus to derive the collision 
terms, we apply the many-body Keldysh-Green function technique 36 combined with the path integral approach. The 
method described below automatically takes into account the many-body effects, namely the Pauli exclusion principle, 
leading to a non-linear form of the collision terms and eventually to suppression of the decoherence. 

The starting point of the derivation is Eqs. (15), (22) for the propagator, in which the integration over the fast 
fermionic fields, "4>l,r, and phase, <j>, has been performed while integration over the phonons and Andreev states 
remains, 



U = J V{X< i }V 2 ne iS/h . 



(78) 



The time evolution in the action follows now along the Keldysh contour 37 Ck, S = J Ck dtL, which goes from — oo to 

+oo, and then backwards. 20,38,39 The interaction is supposed to be switched on and off adiabatically at the remote 
past, t = — oo, and the phonon bath is supposed to be in the thermal equilibrium. The Lagrangian L has the form, 

L = fj {ihd t -Ela z )'n-^2'r*Xqfj(7 x Ti + L ph , (79) 

q 

where L p h is given by Eq. (20). To reduce the time integration along the Keldysh contour to an ordinary time integral, 
we distinguish forward and backward branches of the contour by labeling them with index s = 1,2, and introduce 
the two-component fields 77 s and X*. Since the action is local in time, it can be rewritten as S — dt (L 1 — L 2 ), 
where L s = L[fj s ,rj s , X^]. 

The first step of the derivation is to integrate out the phonon fields, which will give rise to an effective self- interaction 
for the field 77, 

Sint[f)M =~J dtdt' [fj(t)vMt)] S [r z D(t - t')r z ] ss ' [fj(t')a xV (t')} s ' , (80) 
with the kernel D(t — t') given by 

D(t-t') = J2K\ 2 D«(t-t'), D™'(t-t') = -(i/h)(Tc(x^t)X^t'))), (81) 
q 

where Dq(t — t') is the equilibrium Keldysh Green function of the phonons represented with the 2x2 matrix in the 
Keldysh space. In these equations, Tc is the time-ordering operator on the Keldysh contour, and t z is the Pauli 
matrix operating in the Keldysh space; summation over repeated indices is implied. 

The next step is to take advantage of the weak electron-phonon interaction, and to decouple the four-fermionic 
interaction term in Eq. (80) by introducing the Hubbard-Stratonovich field £/^g(i, if), which is a matrix in the 
Keldysh-Nambu-time space. Before doing this, it is convenient to explicitly extract the small parameter \ p h, which 
determines the electron-phonon coupling strength, from the kernel D in Eq. (80) by redefining the kernel, D —> \ p hD 
(^ P h ~ {v 7 ' s){u/bj D ) 2 <C 1, ujd is the Debye frequency). As the result, we get, 

U = JvgV 2 7 1 exp (i/K) j dtdt' ^?7 s (t) [£(t,i')^ -^S(t,t')r z ] ss ' V' \t')+iW[g] , (82) 

where C{t, t') = (ihdt — E*a z ) S(t — t') is diagonal in the Keldysh space, and 

W W=2^J dtdt'G^^t') [T z b(t-t')T z ] ss ' [a x Q(t',t)a x ]y a , (83) 
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Ki(t,t')=ihD ss '(t-t') [a x g(t,t')a x ] S *p . (84) 

Equation (82) describes dynamics of the field r\ interacting with the effective field Q. In terms of the Keldysh-Green 
function for the field 77, 

G$(t,f) = -W)(T C (»£(*)»?}*' (f))>. (85) 

this evolution is described by the Dyson equation, 

(Ct z - T z t[g]r z ) G = I , (86) 

where the self-energy S depends on the effective field Q. A closed equation for Q can be derived by integrating out 
the field r\ in Eq. (82). This procedure leads to the equation, 



u = jvge lS[S]/h , s[gyn = -iTv\n{CT Z -T Z t[g]T Z ) + w[g}, 



(87) 



where Tr denotes both the matrix trace in the Keldysh-Nambu space and the integration over the time variables. 
Noticing that the action is large, S[Q] <~ X~ h , we evaluate the integral in Eq. (87) within the saddle-point approx- 
imation (cf. Rcf. 38). The corresponding saddle-point equation is derived by varying the action with respect to g, 
which yields, 

(Ct z - r z t[g]r z ) g = x ph i. (88) 

Comparing Eqs. (88) and (86), we obtain the relation, g = X p hG. Written in the terms of G, 

(jCr z -r z t[X ph G]T z )G=i, (89) 

the saddle-point equation is the Dyson equation for the qubit Keldysh-Green function, Eq. (85), in which the self- 
energy contains only an undressed vertex part and a free phonon Green function. Including the parameter X p h back 
into the kernel, D, X p hD — > D, we arrive at the expression for the self-energy, Eq. (84), with Q being replaced by G, 
while D is given by Eq. (81). 

To proceed with the derivation of kinetic equation, it is convenient to introduce a triangular form for the Keldysh- 
Green function by performing transformation in the Keldysh space, 

G^LrM-^^ f A ) , L=i=(; "J), (90) 

where G R ^ = G 11 - G 12 < 21 ) is the retarded (advanced) Green function, and G K = G 11 + G 22 = G 12 + G 21 is the 
Keldysh component. Similar relations also hold for the self-energy. Then Eq. (89) takes the form, 

CG = l + tG. (91) 

Kinetic equation is obtained by considering the difference between Eq. (91) and its Hermitian conjugate for the 
Keldysh component, 36 

ih(d t +d t ,)G K (t,t')-E* a [a Zl G K {t 1 t')} = (Yi R G K — G k yA + Y, K G A — G R Y* K ^ (i, t') . (92) 

The right hand side of Eq. (92) describes the qubit decoherence as well as dynamic corrections due to the phonons. 40 
In the absence of the coupling to phonons, the solution of Eq. (92) has the form, 

Gg(t,t') = {-i/h)e- lE >^ t/n F e lE >^' /h , (93) 

where F is a time-independent matrix determined by the initial state of the qubit. When a weak interaction with the 
phonons is switched on, an asymptotical solution to the equation (92) can be written on the form, 

G K (t,t') = (~i/h)e- lE >^/ R \F{{t + t')/2) + F(t,t')] e lE >^' /R , (94) 
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where the matrix F is a slowly evolving on the time scale, %/E* a , function of the global time (t + t')/2, and F is a 
rapidly oscillating small-amplitude correction (see Appendix B). Equations for the matrix elements of F are derived 
in Appendix B, Eqs. (B14), (B15), and for the diagonal matrix elements they read, 

d t F x = -d t F 2 = ~[(2N + 1) (Fx - F 2 ) + F ± F 2 - 1 ] , (95) 

while for the off-diagonal matrix element, F\ 2 = the equation has the form, 

v 



d t F 12 = 



2 (2N + 1 - F z ) + 2i(6 + S Q F Z ) F 12 , F z (t) = (Fi — F 2 ) /2 , (96) 



where v is the phonon-induced transition rate between the qubit levels, 

v = (2n/h) | 7 ;| 2 S(2E: hQ q ) . (97) 

iV = (e 2/3B » — l) 1 is the phonon distribution function at the frequency equal to the qubit level spacing, [3 — 1/kT, 
and S and 5q are small dynamic corrections defined in Eq. (B16). 

For equal times, t = t', G K (t,t) is related to the qubit density matrix, Eq. (40), as follows: G^ a ,(t,t) = 
(— i/Ti) [2p aa >(t) — 8 aa i}, and therefore Eqs. (95), (96) in fact give kinetic equation for the qubit density matrix 
in the interaction picture, p = e' cr * f W i pe~ Mr *- E «*/' 1 , F a — 2p aa — 1, F\ 2 — 2p\ 2 . It is instructive to write Eq. (95) in 
terms of the qubit occupation numbers n a = 1 — p oa , 

d t n 1 = -d t n 2 = -v[(N+l)n 1 (l-n 2 ) - Nn 2 {l-n 1 )\ . (98) 

The right hand side of this equation has the standard form of the electron-phonon collision term, yielding the Fermi 
distribution for the equilibrium occupation numbers, 

This conclusion is consistent with the well known fact that the Fermi distribution of the Andrecv levels gives correct 
magnitude for the equilibrium Joscphson current 28 ' 29 . Furthermore, it follows from Eqs. (98) and (99) that Trp(i) = 1. 
Then equations for the two independent components of the density matrix, p z = (pn — p 22 )/2, and p\ 2 , omitting the 
dynamic corrections, read, 

tPz - -v [ (2N + 1) Pz - p 2 z - 1/4 ] , (100) 

dtPi2 = —v(N + 1/2 — p z ) pi2 . (101) 

These nonlinear equations are drastically different from the linear Bloch-Redfield equation describing decoherence of 
the macroscopic superconducting qubits, 22 and they have qualitatively different solutions, as illustrated in Fig. 6. 
Exact solutions for Eqs. (100), (101) read, 

r (A = Sp z (0)e- rt v 

PA ) l + ^(0)sinh(/?£*)(l- e - r *) ' sinh(/3£*)' [W > 



Pi 2 (0)e-"/ 2 
Sp z (0)sinh(l3E*){l-e- 



PuW = 1 I 7Z TTWZnTT^BZTT, — . ( 103 ) 



where 5p z (t) = p° z — p z (t) is the deviation from the equilibrium, p° z = (1/2) tanh(/3£'*/2). The evolutions of the 
diagonal (relaxation) and off-diagonal (dephasing) parts of the density matrix are qualitatively similar. One may 
distinguish the linear regime, 5p z (0) sinh(/3i?*) <C 1, when the decoherence is determined by the exponential law, 

Sp z (t) - Sp z (0) e~ rt , p 12 (t) = p 12 (0) e- Tt ' 2 . (104) 

However, the decoherence rate, T, becomes exponentially small at temperature smaller than the qubit level spacing, 
(3E* ^> 1. At this temperature, the most interesting is the opposite, nonlinear regime, Sp z (0) smh([3E*) ^> 1. In this 
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case, there is a wide time interval, t <C smh((3E*)/v, where both the relaxation and dephasing follow the power law 
(see Fig. 6), 



£ (>\ 1 - ,.N g!2(0) 1 

Sp,W = - t Pl2 (t) = ^- M -, 



(105) 



and only at very large times, t 3> sinh(/?.E*)/i/, the decoherence undergoes a crossover to an exponential regime 
similar to Eq. (104). We note that the exponentially small relaxation rate in the linear regime is well known for the 
quasiparticle recombination in bulk superconductors. 41 




FIG. 6. Decay with time of a "cat" state (8p z (0) = pi2(0) = 1/2). Bold line indicates evolution of both (normalized) 
density matrix elements for ALQ for 1//3 = Q.2E*; for comparison, exponential relaxation and dephasing of a macroscopic 
superconducting qubit is illustrated with the dashed and dashed-dotted lines, respectively. 

C. Evaluation of transition rate 

We conclude our study with the evaluation of the phonon- induced transition rate, v in Eq. (97). To evaluate 
the transition rate, one needs to specify geometry of the junction in more detail. Let us suppose that our adiabatic 
constriction, Eq. (4), is formed by a hard- wall potential and has an axial symmetry. Under these assumptions, the 
Fourier component of transverse wave function in Eq. (75) has the form, 



F(d±,x) = 2 



Ji(r±(x)q±) 
r±(x)q± 



(106) 



where r± (x) is the radius of the constriction cross section. The magnitude of the relaxation rate essentially depends on 
the parameter r±(0)Q, where r±(0) is the radius in the neck of the constriction, and Q = 2E*/hs is the wave vector of 
phonons responsible for the interlevel transitions; for atomic-size constrictions, this parameter is small, r±(0)Q <C 1. 
Let us assume that the qubit level spacing is not too small, E*/A 3> s/v, then the phonon wave vector Q is large 
compared to the inverse penetration length of the Andreev level wave function, Q ^> k. 

Let us for a moment assume that the Andreev level wave function does not spread out in the electrodes, but remains 
confined in the transverse direction, rj_(x) = const = rj_(0) (see Fig. 7); then the Fourier component in Eq. (106) 
is close to unity, and the interaction region in Eq. (75) is limited by the penetration length of the Andreev state, 
x <~ 1/k, restricting relevant phonon longitudinal wave vectors to small values, q x ~ k <C Q- The transition rate in 
this case reads, 



v = 



7 2 ir 



nQ^R*(Ce/E:)r^(E:) 



(107) 



where t h (E*) ~ E* 3 /hQ 2 D is a bulk electron-phonon relaxation rate at the Andreev level energy (9d is the Debye 
temperature). This result has been derived in Ref. 32 (although neglecting the renormalization effect), and it can 
be qualitatively applied to long constrictions, whose length exceeds the coherence length. For short constrictions 
considered here, L <C l//c, the effect of spreading out of the Andreev level wave function is essential, and the 
approximation, r± = const, is not appropriate. 
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FIG. 7. Interaction region of the Andreev levels with phonons in short QPC (dark shadow) and in long QPC (light shadow); 
in long QPC, increase of the constriction radius (dashed line) can be neglected. 



Let us adopt the following model for the constriction shape, r±(x) = r±(0) (1 — \x\/ L)~ a , a > 1. It is easy to see 
then that the integral in Eq. (75) will have a cut off at \x\ = L <C 1/n. Function J 1 (z a )/z a in Eq. (106) can with 
good accuracy be approximated with the function, (1/2) 9(z — z), where z ss 1, and the integral in Eq. (75) is easily 
evaluated, giving, 



kLvq . 



(108) 



i.e., the transition rate in short constrictions significantly reduces, 
available for the Andreev level-phonon interaction in short QPC. 



This is the effect due to small spatial region 



VII. CONCLUSION 



Let us summarize the outlined theory for the Andreev level qubit (ALQ). ALQ belongs to the family of super- 
conducting flux qubits, but it differs from the macroscopic flux qubits 3-6 in several important respects. First, the 
quantum hybridization of the flux (and persistent current) states in the ALQ loop is produced by electronic back 
scattering in the quantum point contact (QPC) rather than by charge fluctuations on the junction capacitors in the 
case of the macroscopic qubits. Thus, in principle, neither small junction capacitance nor large loop inductance is 
critical for the ALQ operation. Second, the ALQ is based on a QPC with large, almost full, transparency, in contrast 
to classical tunnel junctions employed in the macroscopic superconducting qubits. Large contact transparency is 
required for placing Andreev levels deep within the superconducting energy gap to achieve good decoupling from the 
continuum electronic states. To guarantee good separation of the qubit levels from the continuum, the amplitude of 
the phase fluctuations around the biasing point, <p ~ ^ must be restricted to the small values, which implies small 
inductance of the qubit loop. 

In the tunnel junctions of the macroscopic qubits, Andreev levels are fast variables whose effect, after averaging, 
reduces to the Josephson potential energy added to the Hamiltonian of the loop oscillator. In the transparent junction 
of the ALQ, Andreev levels are slow variables which cannot be averaged out, and the full description includes the 
Andreev two-level Hamiltonian strongly coupled to the quantum loop oscillator. Derivation of the effective two-level 
Hamiltonian goes beyond the tunnel model approximation, which is done by incorporating exact boundary condition 
into the action of the contact. 

The qubit read out is achieved via measuring fluctuating persistent current or induced flux in the qubit loop. To 
simplify the read out, the loop plasma frequency is chosen to be large compared to the qubit frequency so that the 
loop oscillator is "enslaved" by the Andreev levels. Other regimes, e.g. resonance between the Andreev levels and loop 
oscillator, could be also considered, however they remained outside the scope of the paper. Typical circuit parameters 
for the ALQ could be chosen as follows: Ej w A for the single open mode, R < 0.01 with I c <~ 400nA for Nb; 
L -O.lnH, and C ~ O.lpF, giving w p ~ 10 11 sec" 1 , and E L /Ti ~ lO^sec" 1 - 10A Nb /h, while E*/K ~ lO^sec" 1 . 

Full description of the ALQ dynamics and qubit-qubit coupling, including qubit manipulation and read out is given 
by the single-particle density matrix. 

Even if a good decoupling of the qubit states from the continuum electronic states is provided, there are still soft 
microscopic modes in the junction, which could couple to the Andreev levels. These modes present potential source 
of "intrinsic" decoherence, in addition to the commonly considered external decoherence, e.g., due to fluctuating 
biasing and read out circuits. We considered such an intrinsic decoherence of the ALQ related to acoustic phonons. 
It turns out that the collision terms in the kinetic equation are nonlinear, in contrast to the linear master equation 
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for the macroscopic superconducting qubits. 22 This reflects the fermionic nature of the Andreev states and leads 
to considerable enhancement of the decoherence time at small temperature. One can understand this effect in the 
following way. Andreev levels belong to a many-body system of superconducting electrons. Although the macroscopic 
behaviour of the ALQ can be expressed in terms of the single-particle density matrix, the microscopic interaction 
with phonons involves two-particle correlation functions, which are sensitive to the fermionic nature of the Andreev 
states and obey the Pauli exclusion principle. This leads to the reduced probability of the phonon-induced interlevel 
transitions and hence to a slower decoherence. 

Furthermore, the rate of the phonon-induced transitions between the Andreev levels is significantly reduced com- 
pared to the bulk transition rate. The reason is that both the Andreev levels belong to the same normal electronic 
mode; this together with a rapid spreading out of the Andreev level wave function in the contact electrode strongly 
reduces the relevant phonon phase space. 

In the paper, only the case of a single-mode QPC was considered for clarity. However, the approach might be also 
relevant for macroscopic Josephson qubits with tunnel junctions. In junctions with disordered tunnel barriers, the 
open conducting modes with large transmissivity are present. 42,43 This introduces low-energy Andreev levels, which 
implies that quantum phase fluctuations become coupled to these Andreev levels, and the system must be described 
with the effective ALQ-type Hamiltonian. Another kind of qubits where the effective ALQ Hamiltonian might be 
appropriate are d-wave qubits: for the most interesting junction geometries, low-energy Andreev levels, midgap states, 
build up in the junction, which also makes the effective ALQ Hamiltonian to be relevant. 
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APPENDIX A: TUNNEL LIMIT FOR QPC 

In this Appendix, we consider a low transparency QPC, Del, and connect our method to the results of Refs. 18, 19 
for tunnel Josephson junctions. In this limit, the Andreev levels lie very close to the edge of the superconducting gap, 
E a w A, and therefore the field r\ is a fast variable, which has to be integrated out along with all the other electronic 
degrees of freedom. This will result in the effective action for the phase difference alone. After the integration over 
electronic fields, the propagator in Eq. (22) will take the form, 



[4>] = J V 2 r,V 2 ^ L ^e^' h = exp{iS[4>]/h} , (Al) 



S[<t>]/n = -*Tr In (l + + ° g- 1 e^l 2 g e^^j . (A2) 



(1 + Vi?) 2 

Here g is the Green function defined in Eq. (32), 

f+OO 



m = -d/M / °° P e-** - *" + Ag - , (A3) 

J -°° Z7T v /A 2 -n 2 (w + «sgnwO) 2 

and the matrix product in Eq. (A2) includes also the time convolutions. Taking advantage of the small D, and 
expanding the action (A2), the lowest order term reads: 

/+oo 
dhdt 2 g-Hh - h) e-^^/ 2 g(t 2 - h) e *'*(*0/2 , (A4) 
-oo 

where the trace refers to the Nambu space. After taking the trace, the action can be written on the following form, 



/+oo 
-oo 



dtidt 



a(ti - t 2 ) cos 1- p(ti - t 2 ) cos 



(A5) 
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with the kernels a and (3 given by 



a(t) = i{D 12) {A/2hf iH^itA/h) 



(A6) 



I3{t) = -i(D/2) (A/2^) 2 H^\tA/K) 



(A7) 



where floi are the Hankel functions of the first kind. Analytical continuation to the imaginary time in Eqs. (A6) 

and (A7) using relation, K v (t) = (tt /2)i 1+1J (it), leads to the same expressions for the kernels, a and (3, as derived 
in Refs. 18, 19, 34, for tunnel Josephson junctions [with normal resistance of the tunnel junction being replaced with 
the normal resistance of the single-channel QPC: Rn = (-De 2 /wh) -1 ]. 

For a slow varying phase, <f>, on the time scale of the kernel variations, %/ A, the both cosine terms in Eq. (A5) can 
be expanded over the relative time coordinate, r = t\ — t<i, up to the second order, and the effective action takes a 
simpler form, 



S[4>] 



dt 



2 \ 2e 



2e 



L cos ( 



(A8) 



where I c = 2e f_™ dr /3(t) = eDA/2h is the critical Josephson current of the single-channel tunnel point contact, 
and 



5C{cj>) 



- 2 r +oc 3 Dp 1 

J drr 2 [a(r) - (3(r) cos <f>] = __[!- (1/3) cos 0] 



3 De 2 



ft J_ x ^ 32 A 

is the correction to the contact capacitance due to the quasiparticle tunneling. 34 ' 19 



(A9) 



APPENDIX B: DERIVATION OF KINETIC EQUATION 



In this Appendix we derive Eqs. (95), (96) from the equation for the Keldysh function, Eq. (92). It is convenient 
to write the equation in the mixed representation, 

ihd t G* (t) - E* a [a Zl G« (i)] = ( G K - G K E A + G A - G R £ K )^ (t) , (Bl) 

where 

G u {t) = J dr e WT G(t + r/2, t - r/2) , (B2) 
and the products include the convolutions defined by the equation, 

{AB) u (t)=J J dT 1 dT 2 e^-^ + ^-^A^ 1 (t + Tj2)B OJ2 (t-T 2 /2). (B3) 

Neglecting effect of the phonons on the retarded and advanced Green functions of the qubit, G R ' A , we replace them in 
the right hand side of Eq. (Bl) by the free Green functions, G R,A — (fru) — E* a a z ± iO)^ 1 , which are time-independent 
in the mixed representation. Thus, the time-dependence in the self-energy comes only from the Keldysh component, 
G K . The self-energy components take the form, 

<J X , (B4) 



Y$ A \t) = (i%l2)o x 



2^ 



K n R(A) 



D£G 



-n 



D 



R(A) r ,K 



n 



G«-n(t) 



(t) = (in/2) a x J^-(D*-D A )[ (2N n + 1) G*_ n (t) + G*_ n G A _ n } a x , (B5) 

where 
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D 



RA _ f dfl' V iV 



2tt% - fl' ± iO 



, D« = -{i/n)(2N n + l)V n 



(B6) 



Nq = (e^ ftn — l) is the equilibrium phonon distribution function, and T>q is the spectral weight function of the 
phonon bath, 



Pn = 27rsgn(0)^ | 7 ;| 2 S(\Q\ - fi,) . 



(B7) 



After integrating Eq. (Bl) over u>, we obtain equation for the Keldysh function at coinciding times, 

G K (t)=J ^G«(t), (B8) 



ihd t G K (t)~E* a [a Zl G K {t)} =h+h{t)+I 2 {t). 
where Iq = (1/2H 2 ) T> 2E » / h a z , and 
dflduj r+co 



(B9) 



*i(*) = «r 



2h J (2tt) 2 Jo 



dr (2JV n + l)D n ^("-^W^ (t - r/2) - e <(«+n-EW0T a x G^ (t - r/2)a x - /i.e. 



/" dtldw 1 dw 2 r +ca 
J 2(2^)3 X ' 



e i(a»- W1 -n)r (t _ t/2)(Tx G K {t _ t/2) + ^ 



(BIO) 



(Bll) 



For a weak electron-phonon interaction, the right hand side of Eq. (B9) is a small perturbation, which allows one 
to construct an asymptotic solution by using an improved perturbation expansion, 



G% (t) = -(2m/h) e - iE >^ h F u (Xt) + J2 A n *£ n) (*) ^ E> " 



t/h 



(B12) 



71=1 



F u (t)=6{u>-K<'z/h)( F f ) F2 ° (t) ) 



F 12 (t) 
Ff 2 (i) 



(B13) 



In Eq. (B12), A is a formal perturbation parameter, which reflects weak electron-phonon interaction, ~ 2?n, and 
allows one to develop a systematic perturbative expansion. In the zero-order approximation with respect to A, the 
(time-independent) matrix F in Eq. (B13) is the solution of Eq. (B9) without the right hand side. The first-order 
equation determines the time dependence in this matrix (which absorbs formally diverging with time terms in a 
straightforward perturbative expansion), 



OtFr = -d t F 2 = -- [ (2N + 1) (Fi - F 2 ) + F ± F 2 - 1 ] 



(B14) 



d t F 12 = - 



(2N + 1 - F z ) + 2i(5 + S Q F Z ) F 12 , F z (t) = {F 1 - F 2 ) /2 



(B15) 



where v = V 2E */ h /h 2 is the phonon-induced transition rate between the qubit levels, N — N 2E */ n , and quantities 



(1/2^)^(2^ + 1) - 



2ir 



n + 2E*/h 



5 a = (l/2tf)j a[) 



2ir fl + 2E*/h 



(B16) 



determine the phonon-induced shift of the qubit frequency. The higher order equations determine the rapidly oscil- 
lating terms, F^\ in Eq. (B12); e.g., equation for F^ reads, 



(B17) 
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where a± = (1/2) (a x ± ia y ), 

l z = e- iiE ^' h (i//4 - 1<5 ) P? 2 (t) + c.c. , (B18) 



j j* g— 4iE*t/H 



u (N + 1/2) + 2i (6 + S F z )]F 12 (t), S = {l/2% 2 )j^- ^ . (B19) 



It follows from these equations that indeed rapidly oscillates, with the frequency, 4E*/h, and has relatively small 
amplitude, proportional to Kv/E*. 
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